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Abstract 

Shell structure of the single-particle spectrum for reflection-asymmetric deformed cavity is investigated. Remarkable 
shell structure emerges for certain combinations of quadrupole and octupole deformations. Semiclassical periodic- 
orbit analysis indicates that bifurcation of equatorial orbits plays an important role in the formation of this new shell 
structure. 



Theoretical and experimental exploration of 
reflection-asymmetric deformed shapes is one of the 
current topics of interest both in nuclear-structure 
and micro-cluster physics flLpLplH] . 

In theoretical calculations, various approaches like 
Hartree-Fock-Bogoliubov methods, microscopic- 
macroscopic methods and semi-classical methods 
have been exploited for this aim (see Ref. jj] for 
a review). Each method possesses merits and de- 
merits, so that it would be desirable to explore the 
subject with various approaches. 

A basic motive of the semi-classical periodic-orbit 
approach [^[g,0| is to understand the origin of shell- 
structure formation that play decisive role in bring- 
ing about symmetry-breaking in the average poten- 
tial of finite quantum systems like nuclei and micro- 
cluster. If we could understand the origin and obtain 
a global perspective, it would become possible to 
qualitatively predict where we can expect a particu- 
lar deformation to appear in the multi-dimensional 
space spanned by various deformation-parameters 
and the number of constituents of the system. 



system becomes non-integrable when the octupole 
deformation is added, and because the degeneracy of 
the periodic-orbits is then reduced. Contrary to this 
expectation, a significant shell structure was found 
in Refs. ||[| to emerge for certain combinations 
of quadrupole and octupole deformations in the 
reflection-asymmetric deformed oscillator model. It 
was pointed out that this shell-structure enhance- 
ment is associated with bifurcation of periodic or- 
bits. 



It would be very important to investigate whether 
such a mechanism of shell-structure formation is 
special to the harmonic-oscillator model or possesses 
more general significance. 

To explore the possibility that significant shell struc- 
ture emerges in the single-particle spectra for non- 
integrable Hamiltonian, we have carried out an 
analysis of single-particle motions in the reflection- 
asymmetric, axially-symmetric deformed cavity by 
parameterizing the surface as 



In the conventional wisdom, shell-structure would 
be weakened if reflection-asymmetric deformation is 
added to the spheroidal shape. This is because the 



R(9) = R t 




m 2 +m 5 



+ 03^30(0) , (1) 
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where a and b are related with the familiar 
quadrupole deformation parameter S (equivalent to 
S osc in Ref. JkJ) by a = ((3 + 6)/ (3 - 2S)) 2/3 and 
b = ((3 - 28)/ (3 + 5)) 1 / 3 . It reduces to spheroid 
(integrable cavity) in the limit that the octupole 
deformation parameter a 3 vanishes. 

We solve the Schrodinger equation under the Dirich- 
let boundary condition and evaluate the shell en- 
ergy by means of the Strutinsky method. To ef- 
fectively obtain a large numbers of eigenvalues as 
a function of deformation parameters, we have ex- 
amined four numerical recipes; the plane- wave de- 
composition (PWD)pi[, the spherical- wave decom- 
position (SWD)fll2||, the boundary integral method 
(BIM) 0,n4,n5[ and the coordinate-transformation 
method (DIAG)|l|Jl||. The DIAG is the most effec- 
tive method for near-spherical shape, but not good 
for strongly deformed shape. In SWD, PWD and 
BIM, eigenvalue problem is converted to a search of 
the zeros of real function, minima of positive func- 
tion and zeros of complex function, respectively, and 
we found SWD is most convenient for the present 
purpose. Thus we mainly use this method, some- 
times cross-checking the results by other methods. 

Figure [j] shows a typical example of the single- 
particle spectrum calculated as a function of 
the octupole-deformation parameter 03 fixing the 
quadrupole-deformation parameter at S = 0.3. It is 
seen that a new shell structure emerges at about 
03 = 0.2. Deformed magic numbers associated with 
this shell structure are 26,42,70,114,172,. . ., taking 
the spin degeneracy factor into account. Note that 
these numbers appear at intermediate places be- 
tween the magic numbers 20,58,92,138,186,. . . of 
the spherical cavity. This indicates that, due to the 
reflection-symmetry breaking of the cavity, strong 
Al = 3 mixing takes place among levels with large 
orbital angular momenta I in spherical major shells. 

Figure [2] shows shell-structure energies evaluated 
with the standard Strutinsky procedure and plotted 
as functions of the particle number N. As expected 
from Fig. [fl, we confirm here that minima develop 
in association with the formation of the new shell 
structure at about 03 = 0.2. 

To understand the physical reason why such a re- 
markable shell structure emerge for a certain com- 
bination of the octupole and quadrupole deforma- 
tions, and to identify the classical periodic orbits 
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Fig. 1. Single-particle-energy spectrum of the deformed 
cavity plotted as a function of the octupole deformation 
parameter 0,3. The quadrupole deformation parameter 
is fixed at 5 = 0.3. The energy is measured in unit of 
H 2 /MRq, M being the mass. 



(a) 5=0.3, fl 3 =0.0 




Fig. 2. Shell structure energies of the deformed cavi- 
ties with 5 — 0.3 and 03 = 0.0(a), 0.2(b), evaluated with 
the conventional Strutinsky method and plotted as func- 
tions of the particle number N. The energy is evaluated 
by putting R = 1.2(2iV) 1/3 fm and Mc 2 = 938MeV, for 
nuclei. 



responsible for this shell structure formation, we an- 
alyze Fourier transform of the quantum spectrum. 

The single-particle equations of motion for the cav- 
ity are invariant with respect to the scaling trans- 
formation (x,p,t) — ► (x, ap, a _1 t). The action in- 
tegral S-, for the periodic orbit 7 corresponds to the 
length L 7 of it, 

Sry(E = p 2 /2M) = Ipdq = pL j7 (2) 



and the Gutzwillcr trace formula is written as 



p{E) ~ p{E) + J2 A^k^-W 2 cos(fc£ 7 - iTfi-y/2), 



(3) 



where p(E) denotes the contributions of orbits of 
'zero- length', d 7 the degeneracy and p 1 the Maslov 
phase of the periodic orbit 7. This scaling property 
enables us to make use of the Fourier transformation 
of the level density with respect to the wave number 
k. The Fourier transform F(L) of the level density 
p(E) is written as 



F{L) = I dkk-^ d - 2 ^ 2 e- lkL p(E = h 2 k 2 /2M) 

~F(L) + X)^<J(X-L 7 ). (4) 



which may be regarded as 'length spectrum' ex- 
hibiting peaks at the lengths of individual peri- 
odic orbits. In numerical calculation, the spectrum 
is cut off by Gaussian with cut-off wave number 
k c = 1/AL as 
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Fig. 3. Fourier transforms of the quantum level densities 
for the deformed cavities with S = 0.3 and 03 = 0.0(a), 
0.1(b), 0.2(c). The degeneracy index d = 1 (valid for 
generic periodic orbits) and Gaussian cut-off wave num- 
ber fc c = \/300 are used in (|5|). In each panel, lengths of 
classical periodic orbits in the axis-of-symmetry (equa- 
torial) plane are indicated by short (long) vertical lines. 
The lengths are measured in unit of the radius _Ro • 
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(5) 
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The amplitude A 1 (or A") is proportinal to the 
stability factor 1/ y/\2 — Tr M 7 | (in stationary-phase 



approximation), where M 7 is the monodromy ma- 
trix of orbit 7, and expected to be enhanced in the 
vicinity of the bifurcation point where Tr JW 7 = 2 
(see Ref. §). 

Let us investigate how these peaks change when the 
shape parameters of the cavity is varied. Figure 
shows, as an example, how the pattern of the Fourier 
transform (0) changes as a function of 03, fixing the 
quadrupole deformation parameter at 5 = 0.3. The 
highest peaks at the spheroidal limit (03 = 0) are 
associated with triangular and quadrilateral orbits 





Fig. 4. Periodic orbits in the deformed cavity with 
S = 0.3 and 0,3 ~ 0.16 (at bifurcation). For each periodic 
orbit, the length L and the trace of the monodromy ma- 
trix, TrM, are indicated. Those in the axis-of-symmetry 
plane are displayed in the upper panel and those in the 
equatorial plane in the lower panel. Only linear, tri- 
angular and quadrilateral orbits are displayed. In the 
top-leftmost figure, a sphere tangent to the boundary at 
equatorial plane is indicated by a broken line. 




Fig. 5. Bifurcation line of the equatorial periodic or- 
bits in the quadrupole-octupole deformation parameter 
space. For S — 0.3, bifurcation occurs at 03 ~ 0.16. 

in the axis-of-symmetry plane, whose degeneracies 
are two. It is clearly seen that these peaks decline 
with increasing 03. This is because the octupole de- 
formation breaks the spheroidal symmetry and the 
degeneracy reduces to one corresponding to the ro- 
tation about the symmetry axis. On the other hand, 
we can clearly see that new peaks rise with increas- 
ing 03. These new peaks are found to be associated 
with the periodic orbits in the equatorial plane at 
the center of the larger cluster of the pear-shaped 
cavity. 

The key to understand the reason why such equato- 
rial orbits start to play increasingly important roles 
at finite octupole deformation may lie in the fol- 
lowing point: At certain combinations of S and 0,3, 
the curvature radius of the cluster matches with 
its radius, as illustrated in Fig. [|. This means that 



the phase space structure around these periodic or- 
bits locally becomes that of 'spherical' one. Then 
bifurcation of these equatorial periodic orbits oc- 
curs and new three-dimensional orbits are born. Fig- 
ure o shows the bifurcation lines of this kind in the 
quadrupole-octupole deformation parameter space. 
The bifurcation occurs at 03 ~ 0.16 for the case 
5 = 0.3. In general, remarkable shell structure may 
appear along the bifurcation line and stabilizes the 
reflection-asymmetric deformed shapes. 

In conclusion, we have investigated shell struc- 
ture of the single-particle spectrum in reflection- 
asymmetric deformed cavity. It is found that re- 
markable shell structure emerges for certain com- 
binations of quadrupole and octupole deformations. 
Semiclassical periodic-orbit analysis has been car- 
ried out, and it is found that bifurcation of periodic 
orbits in the equatorial plane plays an important 
role in the formation of this new shell structure. 

A more detailed analysis of the equatorial-orbit bi- 
furcation will be reported elsewhere. 
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